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It is now clear that AGN heat cooling flows, largely by driving winds. The winds 
may contain a relativistic component that generates powerful synchrotron radiation, 
but it is not clear that all winds do so. The spatial and temporal stability of the 
AGN/cooling flow interaction are discussed. CoUimation of the winds probably pro- 
vides spatial stability. Temporal stability may be possible only for black holes with 
masses above a critical value. Both the failure of cooling flows to have adiabatic cores 
and the existence of X-ray cavities confirm the importance of coUimated outflows. I 
quantify the scale of the convective flow that the AGN Hydra would need to drive if it 
balanced radiative inward flow by outward flow parallel to the jets. At least in Virgo 
any such flow must be confined to r < 20 kpc. Hydrodynamical simulations suggest 
that AGN outbursts cannot last longer than ^ 25 Myr. Data for four clusters with 
well studied X-ray cavities suggests that heating associated with cavity formation ap- 
proximately balances radiative cooling. The role of cosmic infall and the mechanism 
of filament formation are briefly touched on. 



1. Introduction 

How important are AGN for the dynamics of cooling 
flows? Are cooling flows in approximate steady states, 
or do they evolve significantly over a Hubble time? The 
debate of these issues is more than a decade old, but 
the recent spectacular increase in the quality of the ob- 
servational data has revived it and stimulated a new 
generation of simulations of AGN/cooling-flow interac- 
tion. Here I review these developments. In Section 2 
I summarize the observational situation. Section 3 dis- 
cusses the physics of AGN/cooling- flow interactions; each 
component affects the other and it is important to as- 
sess whether these interactions are stabilizing or oth- 
erwise. Section 4 gives preliminary conclusions from a 
programme of large numerical simulations of jet-induced 
heating that we are carrying out in Oxford. Section 5 ex- 
amines the balance between heating and cooling from an 
empirical point of view and concludes that, when aver- 
aged over ~ 100 Myr, heating and cooling are in balance. 
Section 6 sums up and briefly discusses cosmic infall and 
filament formation. 

2. The story so far 

Data from the Chandra and XMM-Newton missions 
have convinced pretty much everyone that the inter- 
galactic medium in clusters is not a multi-temperature 
medium: although the quality of the fit to the spectrum 
of an annulus can often be improved by including radia- 
tion from plasma at two temperatures rather than one, 
the temperature of the cooler component never differs 
by more than a factor of a few from that of the higher- 
temperature component, and the emission measure of the 
hotter component always dominates by more than an or- 
der of magnitude (Sanders this meeting). Moreover, the 
imaging data show significant inhomogeneities within an- 
nuli that require two-temperature fits, and it seems likely 
that the second temperature component merely com- 



pensates for temperature gradients within each annu- 
lus. These results falsify the picture of distributed mass 
dropout that dominated cooling-flow studies for nearly 
two decades. 

Within a few kpc of the cluster centre, cool gas is 
detected. Ga s at < 100 K is det e cted in the rotation 
lines of CO (iT.azareff et ahl 119891 iMirabel et a,nn98!ll: 
iReuter et a l."l993': l Edgell200l|) . Gas at a few thousand 
Kelvin is detected in the r otation- vibration transitions 
of H2 l)Donahue et a l. 2000: ' Edge et al.ll2002t) . Gas at a 
few times 10"* ^ K is detected t hrough optic al emission lines 
llLvndFjl197nt K^owie et al. i ll 983: Heckm an et 3 11 119891 
IConseli ce et al."200lV Gas at ~ 10^ K is detected in soft 
X-rays (,Fabian et al. 2003b 'l. The spatial distribution of 
the observed emission lines and soft X-rays is consistent 
with the cold gas being organized into filaments, within 
which the temperature increases continuously from the 
centre outwards. The filaments have complex velocity 
fields and it seems very unlikely that they are in dynam- 
ical equilibrium within the cluster potential. The mass 
of gas at say a tenth of the cluster's virial temperature 
is at least an order of magnitude smaller than the mass 
of such gas that was predicted to accumulate over the 
cluster lifetime by the distributed dropout model. More- 
over, is much more centrally concentrated than it was 
predicted to be. 

In parallel with destroying the concept of distributed 
mass dropout, the new data provide clear evidence 
that active galactic nuclei are pumping energy into 
cooling flows. Specifically, the new data show that 
cavities in the X-ray e r nitting gas, like that discov- 
ered by iBohringer et al.l ^ll993^ in the Perseus c l uster 
are common llFa ^iaiT^ tal. I l2000t iBlantqn et al.l 1200 



McNamara et al.ll200lt iNulsen et alJ >2002: H einz et al 
200^ The low X-ray emissivity within these cavities 



must arise because the pressure within them is domi- 
nated either by very hot thermal plasma or by relativistic 
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particles and their associated magnetic fields, or by both. 
Radio emission from many cavities shows that these ob- 
jects are regions of enhanced energy density from ultra- 
relativistic particles and fields, which establishes a clear 
connection between the cavities and activity by AGN 
within the cluster. 

The title of this conference, 'The Riddle of Cooling 
Flows', reflects a widespread feeling that these observa- 
tional results are unexpected and perplexing. I do not 
share this feeling because I have been arguing for over a 
decade (i) that distributed mass dropout makes no sense 
physically and (ii) that cooling can occur only at the cen- 
tre, where it will generate outbursts by the AGN, so (iii) 
there is a fundamental symbiosis between cooling flows 
and AGN. I say this not merely to brag (though I do) 
but to make the point that far from being puzzling, the 
data from Chandra and Newton-XMM are very much in 
line with expectation if one applies sound physical princi- 
ples to systems in which there is gravitationally trapped 
cooling gas. 

A decade ago my student Gavin Tabor worked out de- 
tailed models of cooling flows in elliptical galaxies rather 
than clusters of galaxies because I felt surer of the physics 
in smaller systems: the timescales are very short in el- 
liptical galaxies and the current impact of system for- 
mation and cosmic infall can be much more securely 
neglected. Unfortunately, the quality of the data for 
elliptical-galaxy coolin g flows is st ill very poor [the ex- 
ception s being Gen A (Kraf t et al.ll2003|) and NGG 4636 
(|Kahn et al. 2002')] but I have no doubt that a similar 
picture applies in ellipticals, and strongly suspect that 
the well known correlations between nuclear black hole 
masses and global properties of bulges involves the dy- 
namics of galax y-sized cooling flows in an essential way 
l|RiTmevll2nn8bl) . 

3. AGN as heaters 

We now know for certain that massive black holes lurk 
at the centres of cooling flows. They will accrete ambient 
gas. The natural estimate of the rate at which a black 
hole accretes from an at mosphere of given pressure P 
and sound speed Cg is the iBondil l)l952f) accretion rate 



M 



(1) 



This accretion rate is extremely sensitive to the sound 
speed of the atmosphere near the black hole. This sen- 
sitivity makes it plausible that the black hole is a ther- 
mostatically controlled heater: when radiative cooling 
lowers Cg, the accretion rate rapidly rises and releases 
accretion energy that reheats the atmosphere. Notice 
that the accretion rate depends on Cs even more sensi- 
tively than the fifth power manifest in ^ because P will 
rise slightly as Cs declines. However, the dependence of 
M on Cs through P is weak, and below I shall ignore it. 

3.1. Stability 

Suppose a cooling flow is in a state in which heating 
balances cooling at each radius and consider the stabil- 
ity of this state in both a spatial and a temporal sense. 




Fig. 1. — Pull curve: the cooling coefficient A(r) for a plasma 
with 0.4 of solar abundances. Dashed curves: curves proportional 
to ^-jjat intersect the full curve at T = 3.5 and 4.5 x 10^ K. 



The question with regard to spatial stability is this. Sup- 
pose the heating weakens for some reason, then the net 
cooling that ensues will make the density profile steeper 
than it was. Will the system recover from this excursion 
when the heating rate later increases? The answer to this 
question is 'no' if the AGN merely returns to its former 
power: that power balanced the unmodified cooling rate, 
so it will be less than the enhanced current cooling rate. 
But if the AGN returns to work with renewed vigour 
after its rest, refreshed perhaps by the enhanced Bondi 
accretion rate, then I think it is plausible that the system 
can recover its old density profile provided it heats the 
ambient medium with jets. For in this case the spatial 
distribution of the injected energy depends on the den- 
sity profile of the cooling gas. Enhanced density at small 
radii leads to the jets disrupting closer to the AGN, and 
thus increases the fraction of the AGN's power that is 
dissipated at small radii. If the increased concentration 
of the jets' energy deposition is large enough, gas near 
the AGN will expand faster than gas further out, and the 
density profile will fiatten to its former shape. Clearly, 
this mechanism needs to be quantified by hydrodynami- 
cal simulations. 

To investigate the temporal stability of heated cooling 
fiows, consider the differential equation 
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This equation describes the temperature of a system that 
is heated at a rate that scales ~ [cf. eq. |^] and 

is cooled at a rate that scales as ^ T~'^A{T), as does a 
plasma that is confined by a constant pressure. In Fig- 
ure ^ the two terms inside the big bracket on the right 
of equation (|2J) are plotted, the dotted curves showing 
the first term for two values of the constant A. Thermal 
equilibrium is possible when the dotted and full curves 
cross; for the equilibrium to be stable the full curve has to 
have the larger slope at the point of intersection. We see 
that for A less than some maximum value, stable equi- 
libria exist at temperatures less than that of the peak in 
the cooling curve at 83 000 K. For any value of A at least 
one point of stable equilibrium exists at a higher temper- 
ature. In the case of the upper dotted curve, a unique 
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high-temperature equilibrium lies at 2.7 x 10^ K. For the 
lower dotted curve there are stable equilibria at 10® K 
and 5 x 10® K. The stability of an equilibrium depends 
on the difference in the slopes of the intersecting curves, 
so the high-temperature equilibria are only weakly sta- 
ble in a linear sense, and, in the case of the lower dotted 
curve, can be destabilized by quite modest variations in 
the parameter A. 

The low- and high-temperature equilibria differ 
strongly in the system's accretion rate and luminosity: 
the low-temperature equilibrium has the larger luminos- 
ity by a factor /^/^, where / is the ratio of the temper- 
atures; in the case of the upper dotted curve in Figure ^ 
this factor is 3.4 x 10®. 

Equation is only a toy equation that gives some in- 
sight into a complex dynamical system that is governed 
by non-linear hydrodynamic equations. The tempera- 
ture that appears in the equation should be interpreted 
as the temperature of the accreting plasma at the point 
inside the black hole's sphere of influence where the Ke- 
pler speed equals the sound speed. This temperature 
will not be significantly lower than that, Tmin, associ- 
ated with the Kepler speed at the edge of the black hole's 
sphere of influence. Only in very small galaxies will Tmin 
be less than 10^ K. Hence real accreting black holes may 
not possess the low-temperature equilibrium that the toy 
equation (0) predicts. What the equation does correctly 
predict for systems with small values of A is their insta- 
bility at all temperatures in the range 0.8 to 5 x 10^ K 
and roughly half of the temperatures from there up to 
2 X 10^ K. This wide-ranging instability may well lead to 
runaway growth of the accretion rate and luminosity. 

The coefficient A that determines what equilibrium 
points exist reflects the relative importance of heating 
and cooling. An increase in the mass of the accreting 
black hole within a given cooling flow will be reflected in 
an increase in the best-fitting value of A, and may move 
A to a value at which only a stable high-temperature 
equilibrium is possible. Hence our toy equation suggests 
that black holes above a critical mass (that depends on 
the parameters of the cooling flow) have a single mod- 
erately stable state of low luminosity, while lower-mass 
black holes have one or more low-luminosity states of 
marginal stability and are otherwise liable to experience 
runaway growth of their luminosity. 

3.2. Output channels from AGN 

We still do not have a clear picture of what happens 
when a given mass of gas falls onto a massive black hole. 
Equation ([[J is a plausible estimate of the rate at which 
gas falls from the centre of a cooling flow into the force- 
field of the central black hole. But how much energy is re- 
leased in consequence, and in what form does it emerge? 

The efficiency e = E/{Mc^) with which accretion by a 
black hole releases energy is controversial. Comparison 
between the space density of massive black holes and the 
integrated luminosit y density of quasars im plies e > 10% 
in luminous quasars l|Yu fc Tremainel2002t) . Controversy 
rages as to whether a similar efRciency applies to accret- 
ing systems that are faint in the optical and UV. It has 
long been recognized that implausibly large luminosities 



are derived for many massive black holes if the rate 1^ 
is m ultiplied by an efficiency a. nvwhere near as high as 
10% ijFabian fc Canizaresl llOS^. One school of thought 
argues that around these black holes, which include Sgr 
A* at the Galactic centre and the black holes at the cen- 
tres of all nearby cooling flows, electrons fail to heat and 
the ions carry the accretion energy over the event horizon 
before the electrons can radiate it. To my mind this pic- 
ture, like distributed mass dropout, is physical nonsense 
l)Binnevll2003a|) . A much more plausible explanation for 
the low luminosities of many black holes is that offered by 
Blandford & Bcgclman (1999): a wind from the surface 
of the accretion disk carries off much of the energy and 
angular momentum that is released by accretion. More- 
over, the wind carries away most of the mass that falls 
within the sphere of influence of the black hole at the rate 
given by The upshot is rather a small rate of mass 
flow over the event horizon, and the release of nearly all 
the associated accretion energy as a sub-relativistic wind. 

We are familiar with synchrotron-emitting jets ema- 
nating from the nuclei of the central galaxies of cool- 
ing flows. It is natural to assume that the Blandford- 
Begelman wind is collimated on parsec scales and to iden- 
tify this with the observed jet. But this association is 
dangerous! The jets from M87 and similar galaxies have 
significant bulk Lorentz factors 7, while a wind from the 
accretion disk is expected to be non-relativistic. Conse- 
quently the observed jets probably arise through vacuum 
breakdown in the ergosphe re of a rotating black hole 
ijBlandford Znaiekl I 1 977^ and are not directly con- 
nected to the subrelativistic wind from the disk. The 
observed jets are conspicuous because shocks in them 
readily generate synchrotron-emitting electrons, but they 
may be ephemeral side-shows. In particular, they may 
flicker on and off on the short timescale associated with 
the ergosphere, and we should expect to see systems in 
which the wind is present, but there is no relativistic jet. 
Sgr A* may be just such a system. The connection be- 
tween the mechanical luminosity of the sub-relativistic 
wind and the synchrotron luminosity of the source may 
be weak or non-existent. Jet luminosities derived from 
observations of relativ istic jets, such as that for M87 by 
iRevnolds et alJ l)l996j) constitute lower limits on the total 
mechanical luminosity of the AGN. 

3.3. Collimation 

We should keep an open mind about the extent to 
which the sub-relativistic wind is collimated. The less 
well collimated it is, the more locally it will couple to the 
thermal plasma of a cooling flow. Tabor first investigated 
the case of local coupling th at would apply if the w ind 
were essentially uncollimated (jTabor &: Binnevll993tl . In 
this case the wind would shock near the black hole to pro- 
duce a hot cavity - a similar c avity would be pr o duced 
in th e situation envisaged by ICiotti fc Ostrikeil l)l997l 
l200ll) . in which the flow is heated by inverse-Compton 
scattering of photons. An entropy inversion is created in 
which specific entropy decreases outwards. The region in 
which this occurs is convectively unstable. As in an early- 
type star, convection cells soon eliminate this gradient, 
and the system settles to approximate hydrostatic equi- 
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librium in a configuration that has a constant-entropy 
core. Tabor & Binney constructed models in which this 
constant entropy core joined at some radius to a steady 
coohng-driven inflow. 

The new data have shown that in coohng flows specific 
entropy is not constant near the centre: in reahty the 
mass M{a) with entropy index a = Pp^^ less than some 
value flts the formula 

Af cx (o- - (To)' (3) 

with gf) a constant and e ~ 1.5 ijKaiser fc BinnevI 
l2008albt ) . This finding suggests that the outflow from the 
AGN that heats a cooling flow is coUimated into jets. A 
jet simultaneously heats the ambient medium at several 
different radii, and thus offers an alternative to convec- 
tion as a means of carrying heat from the AGN out to 
the sites of radiative cooling. The observed cavities in 
the X-ray emitting gas provide direct evidence for coUi- 
mation. 

3.4. FR I or FR II sources? 

The radiation of X-rays makes the cluster gas more 
centrally concentrated, and if unchecked will produce 
a central cooling catastrophe i n a typical cooling flow 
within several hundred Myr. iBinnev fc Taboil ijlOQfTT l 
conjectured that at some stage the central AGN 'winds 
back the clock' and restores the cooling flow to the state 
it was in at an earlier time. A few powerful outburst 
could set the clock back several Gyr each, or many weaker 
outbursts could set the clock back tens of Myr each. 

Fanaroff- Riley (FR) II radio galaxies have powerful jets 
that end in a hot spot on the periphery of the source, 
while the weaker jets of FR I sources break up or bend 
dramatically near the centre of the source. FR II sources 
are thought to have kinetic luminosities ^100 times the 
cooling-flow luminosity and lifetimes ~ 100 Myr, so each 
outburst injects energy comparable to that radiated in 
X-rays over a Hubble time. Consequently, they can be 
responsible for reheating cooling flows only if they rad- 
ically reduce the cooling flow's central density to the 
point at which the central c ooling time becomes al most 
as long as the Hubble time . iRevnolds et alJ l)2002l) and 
iBasson fc AlexandeJ l)2003|) have simulated the impact 
on intracluster gas of FR II sources and shown that most 
of the energy is injected well beyond the cooling radius, 
so FR II sources do not reduce the density within the 
cooling radius sufSciently. 

Eilek (this meeting) notes that the great majority of 
sources at the centres of cooling flows are of the FR I 
type. Hence they have relatively small luminosities and 
jets that become unst a.ble at a sm all fraction of the cool- 
ing radius. In Virgo lOwen et~a l. ( 2000) have mapped 
the diffuse synchrotron emission to unprecedentedly low 
surface-brightness levels, and shown that the emission 
has a sharp edge at ~ 0.4rcooi- This finding suggests 
that not only the current outburst, but all its prede- 
cessors have dissipated their energy within this radius. 
Thus we have a consistent picture in which regular weak 
outbursts by FR I sources replenish the energy radiated 
in the innermost part of the cooling flow. Our group 
in Oxford is concentrating on simulating the impact of 



weak FR I sources l)Omma et al.ll200^ . and below I give 
preliminary conclusions from this work. 

The existence and dynamics of FR II sources implies 
that the thermostat on a massive black hole can fail: 
a rapid drop in the central sound speed Cs, caused, for 
example, by a major accretion event, will provoke a pow- 
erful nuclear outburst. However, precisely the power of 
the jets will cause them to burst out of the core of the 
cooling flow and thus strongly diminish their ability to 
heat and expand the gas that is causing the outburst. 
Thus while cooling flows may be stable to small pertur- 
bations, sufficiently large perturbations may cause them 
to flare up dramatically. 

3.5. Local heating versus entrainment 

One can imaging two rather different routes by which 
an AGN might establish a steady state. In the most 
straightforward route it would inject into each radial shell 
heat that balanced the local cooling. In practice a major 
contributor to the heating of the X-ray emitting plasma 
would be the dissolution of cavities. Another possibility 
is that entrainment of ambient gas would establish an 
outffow along the jet axis that balanced a radiatively 
dr iven inflow e lsewhere. The synchrotron map of Virgo 
bv lOwen et Zl ^200ft) cited above limits the applicability 
of this model in the case of Virgo: since outflow appears 
not to extend beyond r ~ 40 kpc, the regime in which 
dissolution of cavities is unimportant cannot be relevant 
outside r ^ 20 kpc, although it might apply at small 
radii. Heating by dissolving cavities must be important 
at radii of order 30 kpc. 

To understand how cavities dissolve, imagine pouring 
a tanker full of olive oil into the ocean. The oil would 
quickly flow over the ocean surface to form a huge sheet, 
only a few molecules thick. In a similar manner the fluid 
inside a cavity spreads out azimuthally when it reaches 
the level at which the ambient medium has a density sim- 
ilar to itself. What makes the process more complex than 
the spreading of oil on water is the absence of a sharp 
density discontinuity in the ambient medium analogous 
to the water/air interface, and ongoing turbulent mixing 
of hot and cold plasma as the cavity rises and dissolves. 

To investigate the regime of balanced in- and out-flows 
that might apply far inside the cooling radius, it is use- 
ful to calculate a different M profile to one that is fa- 
miliar from the discredited mass-dropout picture. Con- 
sider a body of gas, instantaneously in hydrostatic equi- 
librium, in which the specific entropy density s increases 
outwards. In the absence of heating, the rate at which 
mass moves through specific entropy s is 
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dt 



dm ds 
~d7~dt' 



(4) 



I evaluate the first derivative on the right using JSJ, and 
the second derivative follows from the luminosity per unit 
mass 
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Fig-Hlshows the resulting M profile for Hydra. This pro- 
file resembles the M profile of traditional mass-dropout 
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Fig. 2. — The M profile for the Hydra cluster from the data of 
IDavid et al.l 12000) and equations and jFJ. The metallicity is 
assumed to be 0.4Zq. 



theory in increasing with radius, but there are impor- 
tant differences of detail between the two profiles. First, 
interior to r ~ 50kpc the profile of Fig. [21 shows that 



M oc log (r/ro) 



(6) 



with tq ~ 2.5 kpc, while a conventional M profile has 
M (X r. Thus the mass flux rises more slowly with radius 
than in the conventional picture. Second, the profile in 
Fig. 121 reaches values of M that are about twice as large 
as those re ached by the traditional profile for this cluster 
ijDavid et al . 2000). Both results arise because we are 
calculating not the rate at which mass crosses a given 
radial shell, but the rate at which it passes through a 
given specific entropy. Since we have allowed for cooling 
but not heating, the temperature and entropy at a given 
radius are decreasing, so the radius at which a given value 
of the specific entropy is attained is moving outwards at 
the same time that each physical shell of gas is moving 
inwards. Hence, we are determining the mass flux of a 
flow with respect to a moving frame of reference. The 
frame moves fastest near the centre, so that is where M 
is most enhanced with respect to the traditional value. 

If radiatively driven inflow is balanced by outflow along 
the jet, the mass flux carried by the outflow would have 
to increase by ( dAI / ds) ds between the radii at which 
the specific entropy of the inflowing gas moved from s to 
s + ds. If we assume that © holds and that the radial 
density profile can be approximated by a power law (in 
practice p oc r~^'^ in a wide range of r), then one can 
show that 



dM 
~d7 



constant for 3> ctq. 



(7) 



The numerical value of the mass-flow rate along the jet 
axis would be given by Fig. [21 It would be very large 
even at the small radii (r < 30 kpc) at which this picture 
might apply. 

4. The Oxford simulations 

Advances in computer hard- and software finally make 
it feasible to run the simulations of cooling-flow heating 
that I would have liked Tabor to run a decade ago. Then 
we had to settle for spherically symmetric hydrodynam- 
ics, which, as we stressed, inevitably excluded much of 



the k ey physics. We now use ENZO ijBrvan fc NormanI 
Il997l) . a code which does hydrodynamics on adaptive 
Cartesian grids. As far as we are aware, ours is the only 
work on this problem that uses three-dimensional adap- 
tive grids to achieve kpc-scale resolution at the centre of 
a box that encompasses the entire cluster. 

For the reasons given in Section 13.41 we concen- 
trate on low-power outbursts and seek to under- 
stand the dynamics of cavities that are confined to 
within the cooling radius. Previou s simu l ations aimed 
at this problem ("Ch urazov et al.l l200lt iQuilis et alJ 
200lt B riiggen & KaiseJ I200II: iBriiggen et al. I l2002t 
Briiggen fc Kaiseii l200^ have simply deposited energy 



at some arbitrarily chosen location. In our simulations 
ijOmma et al.ll20 0S^ we inject both energy and momen- 
tum in such a way that jets form within the simula- 
tion, and these then heat the ambient medium in two 
dynamically de t ermin ed re gions. In the simula t ions o f 
iRevnolds et^ l)2002() and lBasson fc Alexander! l)2003(l . 
the regions of heating were also dynamically determined. 
However, these simulations differ from ours in two re- 
spects: (i) their jets were > 100 times more luminous 
than ours, and (ii) their jets were imposed through in- 
fiow boundary conditions on a sphere that excluded the 
cluster centre from the computational region. Our algo- 
rithm for jet formation avoids the use of internal bound- 
ary conditions. 

Our main conclusions from this ongoing work are the 
following 

1. The duration of an AGN outburst, rather than its 
power, is what determines the range of radii within 
which its energy is finally deposited. Unless AGN 
outbursts are of rather short duration (< 30Myr), 
cavities move beyond the cooling radius and de- 
posit their energy further from the centre than is 
required if heating is to offset cooling. 

2. Entrainment of cool gas by the jet is an important 
process, as is turbulent mixing of jet-heated and 
ambient gas. 

3. A turbulent vortex that contains a significant quan- 
tity of cool gas trails each cavity. After the jet has 
switched off, the advance of the cavity slows and 
the vortex overtakes it. When it is about twice 
as old as the duration of the outburst, the cavity 
becomes an overdensity. 

4. The overdensity overshoots the radius in the cool- 
ing fiow at which the ambient medium has the same 
specific entropy. Then it falls back and excites the 
strong internal gravity waves in the cooling fiow 
that accompany its azimuthal spread and disap- 
pearance. 

5. Interactions between AGN outbursts are very im- 
portant. If an outburst follows the last sooner than 

70 Myr, the new jet is disrupted at small radii by 
the turbulent wake of its predecessor and no new 
cavity forms. 
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6. Sheets of cold dense gas on the leading edges of cav- 
ities, similar to those observed as X-ray bright rims, 
frequently form as a cavity is driven up through 
cold gas that is rising in the wake of the previous 
cavity. 

5. A statistical steady state? 

The traditional view is that cooling flows are in steady 
states, with cooling balanced by mass-dropout, but a sta- 
tistical steady state could be achieved by a sequence of 
small outbursts by the AGN fjTabor fc Binncv 1993). A 
radically alternative point of view is that cooling flows 
experience a series of cooling catastrophies that provoke 
powerful AGN outbursts, which rapidly inflate the core of 
the cooling flow to a significant extent, setting the scene 
for a p rolonged drift towards the next cooling catas- 
trophe iBinnev fc Tabol Il995t ICiotti fc OstrikeJ Il997l 
l200lt iKaiser fc BinnevI l2003a|l Simulations of the im- 
pact of FR I I sources on cooling flows (Rcvnolds et al] 
120021 feasson fc Alexanded 12003^ imply that powerful 
AGN outbursts of the type we observe locally do not heat 
cluster cores sufficiently to ensure a long interval before 
the next catastrophe. Less strongly collimated outbursts 
may achieve longer intervals between catastrophes, but 
currently a succession of weak outbursts seems the most 
likely possibility. Moreover, recent observations of X-ray 
cavities support this conclusion. 

If we model the energy density within a cavity by that 
of a non-relativistic gas, then the minimum work done by 
an AGN in blowing a cavity of volume F in a medium 
of pressure P is ^PV, of which three fifths resides in 
the cavity. If the fluid within the cavity is relativistic, 
the minimum energy is APV . Since the inflation of cav- 
ities is likely to be highly irreversible, especially in its 
early stages and from the perspective of the ambient 
medium, the actual work done will be larger. Below I 
shall rather conserv atively assume that the work done is 
3PF. iGhiirazov et al. ( 200J^ discuss the physics of cavity 
inflation in some detail. As a cavity is blown, energy is 
carried away into the ambient medium by hydrodynami- 
cal motions, the high-frequency tail of which will be rec- 
ognized as non-linear sound waves. Such wave s may re- 
cently have been observed bv lFabian et al.l (j2003a) in X- 
ray images of the Perseus cluster. The bulk motions asso- 
ciated with lower frequencies may manif est themselves in 
the p eculiar velocities of Ha filaments l)Heckman et al. I 
Much of the energy acquired by ambient gas as 
a cavity inflates will be used to lift gas up in the cluster 
potential. The cooling time of the gas will be incre ased 
but the gas may not actually heat l)Alexandeill2002(l . 

Values for the product PV, the X-ray luminosity of 
the cooling flow Lx (mostly from the classic M value) 
and the characteristic time r = iPV/ Lx are given in Ta- 
ble 1. We see that heating associated with cavities can 
offset radiation from the cooling flow provided outbursts 
occur every 30 to 120 Myr. Simulations show that cav- 
ities rise at of order the sound speed ~ 1000 km s^^, so 
they move through ~ 50 kpc in 50 Myr. Given the high 
frequency with which cavities are found in the cores of 
cooling flows, it is clear that a new pair must be created 



every < 100 Myr. Arguments from synchrotron ageing 
have long implied that more powerful radi o sources have 
lifetimes of this order fredlar et aDll990^ ■ but here the 
conclusion is slightly different: individual AGN outbursts 
might last a significantly shorter time that 100 Myr, but 
the interval between outbursts is of this order. 

The bottom line of this discussion it that, in a sta- 
tistical sense, heat ing by AGN may we ll balanced radia- 
tive cooling, as iTabor fc Binnevl ^^1993^ originally conjec- 
tured. However, the issue cannot be regarded as closed 
because as Burns (this meeting) and Thomas (this meet- 
ing) have pointed out, cosmological simulations imply 
that clusters are accreting relatively cool gas. The power 
required to heat this infalling gas to the virial tempera- 
ture will significantly increase the power an AGN has to 
provide if it is to keep its cooling flow in a steady state. 
The big uncertainty here is how much of the infalling gas 
reaches the clusters radiatively cooling core. Gas that 
ends up outside the core can be largely ignored. 

6. Conclusions 

The Chandra and XMM-Newton missions have opened 
the way a clear understanding of how cooling flows work. 
Distributed mass dropout has been decisively rejected 
and the importance of heating by AGN established. The 
data seem to show that there is approximate balance be- 
tween AGN heating and radiative cooling. Although the 
heating does have to be episodic, the time between heat- 
ing episodes is significantly shorter than even the central 
cooling time, so the radial density profiles of systems do 
not change greatly around a cooling/heating cycle. 

The favoured heating agent is collimated outflow from 
the AGN: a collimated outflow can heat ambient material 
at several radii and values of specific entropy simultane- 
ously. Consequently there is no need for the system to 
develop a constant-entropy convective core. It is proba- 
bly also important for the stability of a cooling flow that 
it is heated by jets: steepening of the ambient density 
profile causes the jets to disrupt at smaller radii. The 
resulting increase in the concentration of the heat source 
may well be successful in counteracting the increase in 
the concentration of the cooling rate that accompanies a 
steepening of the density profile. 

Much of the energy in the heating outflow will be con- 
tained in a sub-relativistic wind. The opening angle of 
this wind may not be small. At its core there may or 
may not be relativistic jets. Observations of synchrotron 
radiation will be largely sensitive to the relativistic jets 
and their fall-out, so such observations may provide a 
very incomplete picture of the overall AGN/cooling flow 
interaction. 

Many questions remain open. Two of the most in- 
teresting, and possibly connected, issues are the roles 
of cosmic infall and cool filaments. Have the fila- 
ments formed through condensation of X-ray emitting 
gas, as has traditionally been argued (|Fabian fc Nulsenl 
Il977|) ? If this conjecture were correct, and in view 
of the thermal stability of gravitationally stratified 
plasma, it is surprising that the filaments seem so far 
from dynamical equilibrium: one would expect con- 
densation to occur onto a rotationally supported disk 
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Table 1 

Parameters for five clusters with cavities. 



System 


erg 


Lx/10*3ergs-l 


T / Myr Reference 


Hydra A 


27 


30 


88 McNamara ct al. C2000'); Nulsen et al. ('2002') 


A2052 


4 


3.2 


122 Blanton et al. (2001) 


Perseus 


8 


27 


29 Fabian et al. (2000): Allen et al. (19921 


A2597 


3.1 


3.8 


79 McNamara et al. ("2001) 


A4059 


22 


18 


119 Huang & Sarazin (1998); Heinz et al. (2002) 



rath er than often radially dire cted filaments. Also, 
as .Sparks. Macchetto fc Golombek ( 1,989.) have stressed, 
one would not expect gas that had condensed from the 
X-ray emitting plasma to have dust that is similar to 
Galactic dust. Perhaps filaments reflect the cosmic infall 
of cool gas and gas-rich galaxies. Very near the centre of 
the cluster this gas might take so long to evaporate that 
it can be lighted up by internal star formation and then 
be observed (Nipoti & Binney in preparation). In this 
case, the AGN would not normally be accreting gas from 
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filaments, but feeding directly on the X-ray emitting gas. 
The AGN's luminosity would then be closely connected 
to the temperature of the coolest X-ray emitting gas, and 
the efficient operation of thermostatic feedback would be 
natural and fairly well modelled by equation 1^)1. 
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